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Abstract 

A class of interacting particle systems on Z, involving instantaneously annihilating or 
coalescing nearest neighbour random walks, are shown to be Pfaffian point processes for 
all deterministic initial conditions. As diffusion limits, explicit Pfaffian kernels are derived 
for a variety of coalescing and annihilating Brownian systems. For Brownian motions on 
R, depending on the initial conditions, the corresponding kernels are closely related to 
the bulk and edge scaling limits of the Pfaffian point process for real eigenvalues for the 
real Ginibre ensemble of random matrices. For Brownian motions on R + with absorbing 
or reflected boundary conditions at zero new interesting Pfaffian kernels appear. We 
illustrate the utility of the Pfaffian structure by determining the extreme statistics of the 
right-most particle for the purely annihilating Brownian motions. 


1 Introduction and statement of key results 


Pfaffian point processes arise in a number of contexts, for example the positions of eigenvalues 
of certain random matrix ensembles (the Gaussian orthogonal ensemble, Gaussian symplectic 
ensemble, the real Ginibre ensemble), and in certain random combinatorial structures such as 
random tilings. In [22] . systems of instantly coalescing, or instantly annihilating, Brownian 
motions, under a maximal entrance law, were shown to be Pfaffian point processes at any 
fixed time t > 0. The aim of this paper is to generalise this result in a number of ways: we 
consider (i) analogous particle systems on Z; (ii) spatially inhomogeneous nearest neighbour 
motion; (iii) general deterministic initial conditions. The Pfaffian point process structure 
survives all of these changes. The key tool in the proof is a Markov duality. 


The introduction is organised as follows: the discrete models are formulated in section 1.1 


the Pfaffian property is stated in section 1.2 the kernels for various diffusion limits are stated 
in section 1.3, applications of the Pfaffian property, and relation to other work are discussed 
in section [T~4l 
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1.1 Mixed models. 


Interaction occurs when one particle jumps onto an already occupied site. In a coalescent 
system there would be an instantaneous coalescence where the two particles merge to leave 
a single particle; in an annihilating system there would be an instantaneous annihilation 
where both particles disappear. We consider the following mixed system, whose dynamics 
are informally described as follows. 

Particle Rules. Between interactions all particles jump independently following a nearest 
neighbour random walk on Z, jumping 

x — y x — 1 at rate q x , and x — 1 —> x at rate p x . 

The parameter 9 E [0,1] is fixed throughout, and when two particles interact they instanta¬ 
neously annihilate with probability 9 or coalesce with probability 1 — 9. 

We denote our particle systems as processes (r] t : t > 0) with values {0, 1} Z , so that rj t (x) = 1 
indicates the presence of a particle at time t at position x. The Markovian dynamics are 
encoded carefully in the generator (see (fl5|)) . We take the left and right jump rates (p x , q x : 
x E Z) to be uniformly bounded, and then this generator determines the law of a unique 
Markov process, for any given initial condition in 77 6 {0,1} Z . We denote its law by P v and 
E, ? on path space with canonical variables ( 77 * : t > 0). 

By choosing 6 = 0 or 6 = 1, our results apply to both purely coalescing and purely annihilating 
systems. We know of three motivations for directly studying the mixed models 6 E (0,1). 


• Excitons in polymer chains. The kinetics of excitons (localised electronic excita¬ 
tions) along polymer chains are observed to exhibit sometimes coalescing collisions and 
sometimes annihilating collisions. Different polymer materials lead to different values 
of the parameter 6 - see Henkel m where three values of 9 are listed for three different 
polymers. 

• Multi-valued voter models. In the multi-valued voter model on Z (also known as 
the Potts model), started from one specific initial condition, the mixed systems arise 
as dual processes. Indeed, starting from product measure, where each of q colours 
has equal 1/q chance, the domain walls that separate regions of sites with the same 
colour behave precisely as the mixed model above with the choice 9 = l/(q — 1). This 
connection is explored in the amazing results of Derrida et al. m , 0 and we discuss 
these further in section [L4l 

• Coalescent mass models. Coalescent mass models are widely studied (see El)- The 
simplest model might be where massive particles perform simple coalescent random 
walks on Z, and where masses add upon coalescence. Consider the case where the 
initial masses are independent at different sites, and chosen uniformly from the set 
{0,1,2,... ,q— 1} for some fixed (integer) q > 1. Then the mass of particles in the 
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future, modulo q, remains uniformly distributed over {0,1, 2,..., q — 1}. Moreover the 
particles whose mass is exactly r modulo q perform a mixed coalescing/annihilating 
system with 8 = 1 jq. This gives only a very partial description, and a full distribution 
of masses at multiple space points remains an interesting problem. 


The spatially inhomogeneous version of the process here occurs in studies on reaction diffusion 
models with quenched disorder (see H3D- The fact that the Pfafhan property extends to 
inhomogeneous nearest neighbour jump rates supports the suspicion discussed in [22] that in 
the continuum the free motion can be any continuous Markov process. 


1.2 Statement of main theorem 


We have decided to recall the definition of a Pfafhan point process on Z. A variable y with 
values in {0,1} Z can be thought of as a point process on Z, with y(x) = 1 corresponding to a 
particle at x. In this context, the definition of a Pfafhan point process (see [19]) is as follows: 
there exists a matrix kernel K : Z 2 —> M 2x2 , so that for any n > 1 


B[ V (x 1 )...r,(x n )] =Pf (K ( Xi , Xj ):i,j < n) for any disjoint aq,..., x n . 

Pf (A) is the Pfafhan of an antisymmetric matrix, and the appendix contains definitions and 
the simple properties of Pfafhans that we use. The matrix kernel K can be written as 


K (y,z) 


Kn (y,z) Ki 2 (y,z) \ 
K 2 i (y,z) K 22{y,z) J 


for K ij : Z 2 —> M, and these entries must satisfy the symmetry conditions 


K = -K ji(z,y) for all i, j E {1,2} and y,z £ Z 

which ensure that the 2 n x 2 n matrix (K(aq, Xj) : i,j < n), built out of the n 2 two-by-two 
blocks, is antisymmetric. 


The matrix kernel for our interacting particle system is constructed from a single scalar 
function : y < z) dehned as follows. For y,z £ Z with y < z, and for rj E {0,1} Z , 

we define 

v[y,z)= r /( x ') if 2 /<A 

y<x<z 

and rj[y, y) = 0, and we dehne the ’spin pair’ by 

OyM = {-0)^ Z) . 

We use the convention that 0° = 1 so that when 0 = 0 the spin reduces to the indicator of 
an empty interval, that is cr y , z {ri) = I(r/[y, z) =0). We now set 

K t (y, z) = E, ; [<J y , z {r}t )\, for t > 0, y, z E Z with y < z. (1) 
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We also need the difference operators D + and D , defined for / : Z —> M by 
D + f(x ) = f(x + 1) - f{x), D~f(x ) = f(x - 1) - fix). 


Theorem 1 For any initial condition r] E {0, 1} Z , and at any fixed time t > 0, the variable 
T) t is a Pfaffian point process on Z with kernel K given, for y < z, by 


K iy,z) 


1 

i + e 


K t (y,z) -D%K t (y,z) 

— DfK t (y,z) D+D+K t iy,z) 


( 2 ) 


and Kio(y,y) = D^Ktiy, y), ond other entries determined by the symmetry conditions. 
(The notation Df means that the operator is applied in the ith variable.) 


Remarks. 


1. Random initial conditions. For random initial conditions the law of ry is not in general 
a Pfaffian point process, though by conditioning on the initial condition rj o the correlation 
functions can always be expressed as the expectation of a Pfaffian with a random kernel K/'?o) 
depending on t]q: 

p ( t n \x 1 ,...,x n ) = E(Pf (K no {x i ,x j ) : i, j < n)). (3) 

For certain random initial conditions, including the natural case when the sites ( rio(x ) : x E Z) 
are independent, the expectation can be taken inside the Pfaffian and the process does remain 
a Pfaffian point process See the first remark at the end of section [2j 


2. Thinning. The parameter 0 enters into the kernel only as a scalar mulitiplier. Instantly 
coalescent systems and instantly annihilating systems are related by a well known thinning 
relation (see [22] section 2.1 or Cl). This extends to our mixed systems as follows: consider a 
two colour system of particles, red R and blue B, that move independently between reactions, 
and at reaction times (when one particle lands on top of another) transform via the rules 


R + R-+R, R + B^B, B + B^ 


R with probability 8, 

B with probability 1 — 9. 


Note that the full system of particles, where one ignores colours, is a coalescing system, but 
the blue particles alone follow the mixed model we study. However, if initially particles are 
coloured 

blue with probability and red with probability 

then this property is preserved at all subsequent times. Indeed one can check that each single 
collision preserves property. Thus at any fixed time the mixed system is a thinning, by the 
factor 1/(1 + 9) of the full coalescing system. Thinning also acts naturally on Pfaffian point 
processes, changing the underlying kernel by the same factor. However this connection seems 
to relate the two systems only when the initial conditions are similarly related by thinning, 
and so does not apply to the deterministic initial conditions stated in the Theorem. 
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2. Strong thinning. The Pfaffian point process with kernel ([ 2 ]) but where the scalar 
multiplier 1/(1 + 9) is replaced by an arbitrary A E (0,1) also arise. When A E (0, |) 
these correspond to thinning a coalescing system by more than the factor needed to reach 
annihilating systems. Consider a two colour system as described above, but with reactions 


R + R^-R, B + B->B, B + R^f 


R with probability 
B with probability ^. 


If we initially choose colours independently as 


blue with probability A and red with probability 1 — A 

then this property is preserved at all later times. Thus the sub-population of blue particles 
is a thinning, by the factor A, of a coalescing system and hence Pfaffian with the coalescing 
kernel thinned by the factor A. 


1.3 Continuous limits 

The scalar function Kf(y , z) that underlies the Pfaffian matrix kernel K can be characterised 
as the solution to a system of differential equations indexed over part of the lattice. We define 
a one-particle generator L, acting on / : Z —> M, by 

Lf{x) = q x D + f(x) +p x D~ f(x). (4) 

The intuition is that L is the generator for a single dual particle. We will show that the 
function ( K t (y , z) : t > 0, y, z E Z, y < z) is the unique bounded solution to the equation 

( d t K t (y , z) = ( L y + L z )K t (y , z) for y < z, t > 0, 

< K t {y,y) = 1 for all y, t > 0, (5) 

{ K 0 (y,z ) = <Ty tZ {rj) for y < z. 

(The notation L y is used to indicate that the operator L acts on the y variable.) 

This differential equation characterisation ([5]) lends itself naturally to asymptotic analysis, 
where its large time and space behaviour is determined by a similar limiting continuum kernel 
I\[ c \y,z) solving a continuum PDE. We can in several natural cases solve these limiting 
continuum equations explicitly, and therefore add to the growing zoo of concrete known 
kernels for Pfaffian point process on M. 

In the cases below we use the diffusive scaling 

Xt e) = Ve-n(x)6 ex ( 6 ) 

and we check, at a fixed t, that —y Xt as e —> 0 (considered as random locally finite 

point measures on M with the topology of vague convergence). By choosing p x ,q x ,9 and the 
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initial conditions we establish the Pfafhan property for various continuum systems at fixed 
times. In each case the limiting point measure Xf is a Pfafhan point process on with a kernel 
of the form 


K {c \y,z) 


1 

1 + 6 


K[ c \y,z) -d z K ( t c \y,z ) \ 
~~d y Kf C \y, z) dy Z Kf C \y, z) )' 


and K^(y,y) = ^ d 2 K t (y,y). 


for y < z, 


( 7 ) 


We record some specially chosen cases where ( y , z) can be found explicitly. The limits can 

also be identihed as the law at time t for a suitable system of reacting Brownian particles. 
In section [3] we will recall the idea of the maximal entrance law for these systems (which 
can be constructed in each case as the limit of Poisson initial conditions as the intensity 
increases to infinity). We emphasise that the Pfafhan property holds for all deterministic 
initial conditions - but the maximal initial condition leads to a simple initial condition of the 
PDE’s determining the kernel and hence simple explicit solution formulae. More importantly, 
as explained in [22] for purely coalescing or annihilating systems, the maximal entrance law is 
distinguished in that the solutions are then invariant under diffusive rescaling. Moreover the 
solutions from a large class of initial distributions becomes attracted under diffusive rescaling 
to the solution started form the the maximal entrance law. 


Theorem 2 Fixt > 0 throughout. Recall erfc{x) = 1—er/(x) for erf(x) 


2 


Jo exp( 


u 2 )du. 


(A) Brownian motions on M - maximal entrance law. 

Take p x = q x = 1 for all x and the initial condition rjo = 1. The diffusion limit Xt has 
the Pfafftan kernel 0 on M with 

Kf>(y,z) = erfc (^|) . 

Moreover, this is the kernel for mixed coalescing/annihilating Brownian motions on M 
at time t under the maximal entrance law. 


(B) Brownian motions on M - half-space maximal entrance law. 

Take p x = q x = 1 for all x and half-space initial conditions: r]o(x) = 1 for x < 0 and 
Vo(x) = 0 for x > 0. The diffusion limit X t has the Pfaffian kernel 0 on M with 

tAc), x , f 27 tf^/i(u — v) ( u-v) 2 fu + v\ , , 

A, (v, *> = !+/ dud,, (9) 

2 vt 

Moreover, this is the kernel for mixed coalescing/annihilating Brownian motions on M 
under the half-space maximal entrance law. 
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(C) Killed Brownian motions on [0,oo) - maximal entrance law. 


Take 


{ 1 for x > 1, 

2 for x = 0, 

0 for x < 1, 


( 1 for x > 1, 
\ 0 for x < 1. 


Also take rjo(x) = 0 for x < 0 and r]o(x) = 1 for x > 1. T/ie diffusion limit Xt has the 
Pfaffian kernel 0 on (0, oo) with 


AfW) = l-er/(^) erf 


z - y \ 

2y/21) 


( 10 ) 


Moreover, this is the kernel for mixed coalescing/annihilating Brownian motions on 
[0, oo), killed at {0}, under the half-space maximal entrance law. 

(D) Reflected motions on [0, oo) - maximal entrance law. 

Take 

_ ( 1 for x > 1, _ f 1 for x > 1, 

\ 0 for x <1, ^ x \ 0 for x < 1. 

Also take rjo(x) = 1 for x > 0 and rjo(x) = 0 for x < 0, The diffusion limit X t has the 
Pfaffian kernel <[?|) on [0, oo) with I\\. c \y,z) given by 

1 - J /r du !l dv i erf " (^tt) erf & r) + er/ " (^tt) er/ (^7r)) - 

2\/t 2\/t 

O') 

Moreover, this is the kernel for mixed coalescing/annihilating Brownian motions on 
[0, oo), reflected at {0}, under the half-space maximal entrance law. 


The final two examples illustrate the fact that we allow spatial inhomogeneities in our discrete 
models. One point intensities pt(y) = —d z K[ c \y, z)\ z=y can then be read off from the Pfaffian 
kernels yielding for the absorbing case 


Pt{y) 


erf (- 


1 + 0^2fft \>/2t 


and for the reflecting case 


Pt(y) 


1 

T+0 




+ 



y 2 

4 1 erfc 



The figure below illustrates these for the coalescing case (9 = 0atf = l/4. Note that both 
intensities converge as y —> oo to that of example (A), the free Brownian motions. 
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Figure 1: Particle density at time t = 1/4 for coalescing Brownian motions killed at the origin 
(example (C)), and reflecting at the origin (example (D)). For the absorbing case, particle 
density vanishes at the origin due to the trap at 0. 

1.4 Related work and applications. 

The real Ginibre random matrix ensemble. The real Ginibre random matrix ensemble 
is the N x N random matrix formed by taking independent N( 0,1) Gaussian variables as 
entries. The positions of the eigenvalues form a Pfafhan point process (see mi, mi). 
Letting IV —> oo there is a limiting kernel describing the eigenvalues in the bulk of the 
spectrum. The paper (22], thanks to the comments of a referee, noted that the kernel for 
the real eigenvalues in the bulk can be written in the form 0 where Kt(y, z) is given by Q 
when 6 = 1, pure annihilation, and t = 1/4. 

A similar large N limit can be taken for the real eigenvalues in the real Ginibre ensemble 
near the right hand edge of the spectrum. Again the limiting kernel can be written in the 
form Q where Kt(y,z) is now given by ([9]) when 6 = 1 and t = 1/4. (The original version 
of [4] has a slightly incorrect derivation of one of the four entries in the limiting edge kernel 
K 77 */e >re {'!J, z )- However this can be easily corrected and all four entries agree with the kernel 
above (see the erratum [5])). We state these observations as a corollary. 


Corollary 3 The large N kernels for the real eigenvalues of the real Ginibre ensemble, in 
the bulk or at the edge, agree with the kernels for annihilating Brownian motions, at time 
t = 1/4, started from the maximal entrance law or the half-space maximal entrance law. 


Gap probabilities. A standard question for the determinantal and Pfafhan point processes 
arising from random matrix models is the problem of gap probabilities 

p(t, [L , R\) = P[there are no particles in an interval \L, i?] at time t] 









and their asymptotics as R—L —> oo. A particularly interesting special case of gap probability 
is p(t, [L, oo)) under the half space initial condition, which gives the law of the right most 
particle. 

These probabilities can be expressed in terms of the Pfaffian kernels as Fredholm Pfaffians, 
and their study for the Pfaffian point processes arising from classical random matrix ensembles 
(GOE and GSE) is described in [2TJ. For the interacting particle systems described here, 
bulk gap probabilities were studied by Derrida and Zeitak [8]. They were motivated by 
the connection with the Potts model described above, and this led them to analyse the gap 
probability under a product Bernoulli initial condition, in the translation invariant case. The 
leading term of their asymptotic might be expressed as 

lim p(t, t l / 2 [L , R\) « e -m(R-L) for R _ L i arge (12) 


where 



Although they do not identify a Pfaffian point process, the arguments in [8] exploit Pfaffians, 
and manipulations similar to those exploited for random matrix ensembles, and this paper was 
a motivation in our search for an underlying Pfaffian point process for these mixed models. In 
a subsequent paper we will show that these techniques, starting from the Fredholm Pfaffian, 
will confirm the asymptotic (12) for all 9 E [0,1] for Bernoulli initial conditions, and then 
extend them to more general initial conditions. 


One non-translation invariant case has already been investigated, namely the edge gap prob¬ 
ability p(t, [L, oo)) for the annihilating particle system. This corresponds to the largest real 
eigenvalue in the real Ginibre ensemble and has been studied in m and m- The main 
theorem in m can be expressed as follows. 


Theorem 4 Consider the system of instantaneously annihilating Brownian motions on M 
with half-space maximal initial conditions. Then 


Pit , [L, oo)) 


1 - er/c(^) + 0(e «) for ± -> oo, 

ex P^^i)A + 0(1)) for 7t^-°°- 


(13) 


The techniques used here start with the Pfaffian point process kernel, but exploit a repre¬ 
sentation for the terms in the Fredholm Pfaffian in terms of a discrete time simple random 
walk, reducing the asymptotic to a simpler asymptotic for a single random walk. We believe 
the same techniques will apply for the mixed coalescing/annihilating models. 

It is interesting to compare the above tail gap problem with the results connecting the statis¬ 
tics for interacting particles in the KPZ universality class, started with half-space initial 
conditions, with the Tracy-Widom distribution corresponding to the behaviour of the largest 
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eigenvalue for Hermitian matrix models. The fact that annihilating Brownian motions ap¬ 
pear as the limit for a large number of discrete interacting particle models, suggests that 
the distribution of the largest real eigenvalue for the real Ginibre ensemble may define an 
interesting universality class associated with statistics for non-Hermitian matrices. 

The underlying Pfaffian structure should be useful in the study of many questions concerning 
these models, for example multi-time correlations. The extended Pfaffian property found for 
annihilating systems in |24j was shown in |il2j to extend to for the discrete space inhomoge¬ 
neous process studied here, solely in the annihilating case. The persistence exponents found 
by Derrida, Bray, and Godreche |7j suggest that Pfaffian structure will be useful for time 
dependent properties as well. 


2 Proof of the main result. 


We start with a summary of the main steps, which follow similar lines as [22]. The key tool is a 
Markov duality. Indeed for any n > 1 the product of n spin pairs y —> a yitV2 ( 77 )... o~ y2ri _ ] , y2n (?/) 
is a suitable Markov duality function, as shown in Lemma [5j Exploiting this allows us to 
calculate the expectations 

\ a yi,y2(vt) ■ ■ ■ (T J/2n-l,J/2n (Vt)] 

as the solutions of 2n dimensional (spatially inhomogeneous) lattice heat equations. This is 
similar to the Markov dualities used in [3] to study the ASEP and q-TASEP models. The 
dual process can be taken to be (a spatially inhomogeneous version of) the one dimensional 
Glauber spin chain (see the Remark after Lemma [5]). This model is known to be solvable by 
mapping to a system of free fermions operators (see Fermions are naturally associated 

to Pfaffians, and it turns out that the duality expectations are given by 2 n x 2n Pfaffians 
of a matrix built from a scalar kernel Kt(yi,yj), as shown in Lemma |6j The final step is to 
reconstruct the particle intensities from the product spin expectations. This is possible via 
the identity 

»(»)=(14) 

This leads to a linear reconstruction formula for the n point intensity in terms of the 2 n x 
2n Pfaffians for the product of n spin pairs. The Pfaffian structure is preserved by the 
reconstruction formula, but the matrix breaks into 2 x 2 blocks corresponding to the spin 
pairs, and this yields the desired matrix kernel K (yi,yj). 

We now present the details. The generator of the process is given, for F : {0,1} Z —> M that 
depend on finitely many coordinates, by 

£F( V ) = ^q40nV a x , x -i) + (l-0) F (v C x , x -i)- F (v)) 

ceEZ 

+ X> {0F(r,Z- l,J + (1 - 0)F{rf x _ 1>x ) - F( V )) , (15) 

x(EZ 
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where yf. y (respectively 'ff x . y ) is the new configuration after a jump from site x to y followed 
by instantaneous annihilation (respectively coalescence). These are defined, when x 7 ^ y, by 

Vx,y( Z ) = vZ,y(z) = v(z) for 2 (jL {x, y}, 
i Vx,y( x ) = r&,y(x) = 0, 

Vx,y(y) = iv(x) + v(y)) rnod(2), 

. Vx,y(y) = min{l, 77 ( 2 )+ 77 ( 2 ;-i)}. 

The key to the argument is the following Markov duality function. For n > 1 and y = 
(yi,, y 2 n) with y\ < yi < ... < yin we define the product spin function by 

n 

=Yl^y2i-i,y2i(v)- 

i =1 

Note that £ y (j 7 ) depends only on finitely many coordinates of 77 and so lies in the domain of 
the generator C. The Markov duality is encoded in the following generator calculation. 


Lemma 5 For yi < yi < ■ ■ ■ < Uin the action of the particle generator L on £^( 77 ) is 

2 n 

2—1 

where L yi , given by 0 - acts on the coordinate yi in T, y . 

Remark. We do not make use of a dual Markov process, but this lemma could be cast into 
the standard framework (see Ethier and Kurtz [9| chapter 4) relating two Markov processes. 
The dual process can be taken to be a finite system of particles with motion generator L 
that are instantly annihilating (with state space the disjoint union U(( l= gM 2m ). However 
this annihilating system describes the motion of domain walls in (a spatially inhomogeneous 
version of) the Glauber dynamics for the Ising spin chain [T3] and the dual process could also 
be taken to be this spin chain. The formulae connecting a set of spins (<r(x) G {—1, +1}, a; G 
Z) to the positions a domain wall separating different spins at x and x + 1 are 

V(x) = -(“l) m ,2/J = cr(x)a(y). 

We do not exploit the link between the spin chain and annihilating systems but it is the 
origin of our use of the term ’spin-pair’. The fact that the dual process both for coalescing or 
annihilating systems can be taken to be the spin chain perhaps explains why both systems 
are Pfaffian. 

Proof. Each term in the generator C involves a modified configuration, 77 “ or 77 £ , which 
differs form 77 on at most two neighbouring sites. The condition that yi < yi+i ensures that 
this modified configuration will agree with 77 on all but at most one of the intervals [yii -\, yn), 
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and hence the value of at most one of the spin pairs a y2i _ 1 ,y 2z will change. This allows us to 
separate the action of the generator as follows 

n / n 

£ s y( ? ?) = ]C n a m-iMjiv) 

i=1 \j= 

We turn our attention to a single spin pair cr ytZ (ri). Fix y < z and consider the part of the 
generator 

Y & ( 9cr yAVx,x-l) + (1 - 9 AyAVx,x-l) - ^yAA) 

corresponding to left jumps. The terms in this sum indexed by x < y and by x > z+1 are zero, 
as the modified configuration is unchanged in the interval [y, z). The terms corresponding to 
y+l<x<z—1 are also zero since we claim that 

6(7 yAVx,x- 1 ) + (1 - 6)<TyAVx, x -l) ~ °yAA = 0 - 

Indeed, since {x — l,x} G [y,z), the left hand side is proportional to 

0 )+V(x) mod(2) _|_ ^ _ Qj (_Q^mm(l,T)(x-l)+r)(x)) _ (_Q^ri(x-l)+T)(x) 

and checking the three cases r/(x — 1) + rj(x) 6 {0,1, 2} shows that this is always zero. (This 
identity is similar to a key quadratic identity (relating the exponential parameter r to the 
ASEP parameters p, q) that lies behind the ASEP dualities in section 4.1 of [3].) Thus jumps 
between sites both lying outside or both lying inside the interval [y, z) give zero contribution 
to the generator and only the terms when x = y or x = z, where a jump crosses an endpoint 
of the interval, contribute. It was this key property that we looked for when trying to find a 
duality function. 


C o 


V2i-l,y2i 


{A- 


(16) 


For the two remaining terms, when x = y or x = z, two similar short calculations lead to 

9°yAVy,y- 1 ) + (1 - Q)°vAVy,y-l) ~ °yAV ) = °y+lAA ~ a V AV ) = Dy°yAA, 
9cr vA r l a z,z- 1 ) + (! - Q)°yAv C z,z- 1 ) - a vAA = <7y,z+i(v) ~ °yAA = D+a y Av)- 

Collecting the terms from all possible x gives 

Y qx ( 6(J yAVx-i,x) + (1 - VAyAVx- I,x) - a yAA) = %D y °yAA + q z D+a y , z (r]). 

Similar calculations for the terms corresponding to right jumps show that 


Y Px ( 6a yAVx-i,x) + (1 - 9 AvAVx-l,x) ~ °vAA) = PyDy ayAA+Pz D z °yAA 


and hence 


£'@y,ziA) — (Cy T Lz) &y,z(A’ 


Using this in (16) completes the proof. 
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For y = (yi, ..., y 2n ) with yi < ... < y 2n , the expectation u(t , y, y) = E r/ [E^ry)] satisfies 


dtu(t, y, rj) = £u(t, y, rj) = E,, [CT, y {y t j] = E r/ 


2 n 

Yi Ly i 'Ey( r lt) 

_i=1 


2 n 

2—1 


Indeed the first equality is the Kolmogorov equation for the Markov process 77 , the second 
equality is due to the Markov semigroup and its generator commuting, the third equality is 
due to Lemma [5} Thus the duality allows us to recast an infinite dimensional Kolmogorov 
equation in ( t , rj) as a finite dimensional ODE in (t, y). We now show that this ODE is exactly 
solved by a Pfaffian. 


Lemma 6 For all y £ {0, 1} Z , for all n > 1, yi < ... < y 2n and t > 0, 

E„ [Xy(r k )] = Pf(K<- 2n '>(t,y)) 

where K^ 2n \t,y) is the anti-symmetric 2 n x 2n matrix with entries K t (yi,yj) for i < j for 
the function Kf defined in 0- that is K t (y,z ) = E v [a y , z (y t )\. 


Proof. For n > 1 denote 


F 2n 

= {y £ Z 2fc : yi < .. 

■ ■ < 2/2n}, 

v 2n 

= {y £ I? k : yi < .. 

■ • < 2/2 n}, 

dV 2 ( S 

= {y £ Z 2n : yi < . 

■■<y%= 2/i+l < 2/7+2 < • • • y 2 n}, 

dV 2n 

— U 2n ~ 1 dV < ' l ' ) 

~ u i= 1 uv 2n ■ 



We now detail a system of ODEs indexed by V 2n , which will involve driving terms indexed 
by dV 2n . Fix an initial condition y £ {0, 1} Z and n > 1 and define 

■u ( 2 n) (i,y) = E^ [T, y (y t )\ for t> 0 , and y £ V 2n . 

For n > 1 , tt( 2n ) solves the following system of ODEs: 


(ODE) 2n 


dtu (2n \t,y) 
< (i, y) 

u ( 2 )i) (Qj y) 


E^i L y .M 2n \t,y) 

u (2n-2)(^ y*,*+l) 


on [0, 00 ) X K 2n , 
on [0, 00 ) x <91^, 
on l/ 2n . 


The notation y ,,l+1 is for the vector y with coordinates y^ and y^+i removed. Thus, when 
n > 2, for y £ we have y v+1 £ V 2n _ 2 . (ODE) 2n is a system of ODEs indexed over V 2n . 

For n > 2, to evaluate L yi u^ 2n ^ one may need the values of u^ 2n l at some points y £ dVf ^, 
which then act as driving functions for the differential equation. The second equation, which 
we call the boundary condition, states that these can be deduced from the values of u^ 2n ~ 2 \ 
Indeed the boundary condition follows simply from the fact that 


E y (y) = Eyi.i+i( 77 ) for y £ dV.^ and n > 2. 
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By setting u ^ = 1, we may suppose the equation holds also for n = 1, encoding the fact 
that u( 2 \t, (y, y)) = 1 for all t > 0, y £ Z. 


The infinite sequence of equations ((ODE) 2n : n = 1,2,...) are uniquely solvable, within the 
class of continuously differentiable functions satisfying sup f>0 sup yg y 2n \u^ 2n \t, y)\ < oo. In¬ 
deed the boundary condition for u ^ is simply that u^ 2 \t, (y, y)) = 1 , and standard (weighted) 
Gronwall estimates show uniqueness of solutions of (ODE) 2 . Inductively, the boundary con¬ 
dition for u ty2n ' > is given by the uniquely determined values of 4 j( 2n = 2 ) and hence u l ' 2n ' ! can be 
found uniquely from (ODE) 2n . 


We now check that (Pf (K^ 2n \t,y)) : n = 1,2,...) also satisfies ((ODE) 2n : n = 1,2,...). 
First we consider the initial conditions. Fix y £ Ejn and choose yo < yi- For 9 > 0 the entries 
in the Pfafhan at time zero can be rewritten as 

Ko(Vi,W)=<W>?)= 

The Pfafhan identity (see appendix) 


shows that 


Pf [ — : 1 < i < j < 2n ) = 1 — - — 1 when a* / 0 for all i, 

0204... 02n 


pf (a-< 2 "> ( o, 9 ))= n ( , Tlrtl 1 = 


(17) 


II (— 9^v[yo,y2i) 

By letting 9 4- 0 we hnd the same hnal identity is true when 0 = 0. 


Next we check the boundary conditions. We fix y £ dV^, t > 0 and write K 1 - 2 ^ for 
K^ 2n \t,y). By conjugating with a suitable elementary matrix E, the matrix 

j^( 2 n) _ E T x( 2n )E 


is the result of subtracting row i + 1 from row i, and column i + 1 from column i. The 
Pfafhan identity Pf ( E T AE) = det(F)Pf (T) ensures that Pf (K <2n ^) = Pf (lf( 2n )). However 
the equality y* = iji+i implies that the zth row of has all zero entries except for 

= 1- Performing a Laplace expansion (see appendix) of the Pfafhan of K < ' 2n '> along 
row i shows that, when n > 2, 

Pf (K^ 2n \t, y)) = Pf (K^ 2n ~ 2 \t, z)) where 2 = y^ +1 . 

When n = 1 we obtain that Pf (K 1 ' 2 ' 1 ) = 1. This is exactly the desired boundary condition. 

Finally we check the differential equation in (ODE) 9n . The entries in K^ 2n ' 1 (t, y) solve (ODE) 2 , 
that is d t K t (y u yj) = (L yi + L Vj )K t (yi,yj). The Pfafhan Pf (K^ 2n \t,y)) is a sum of terms 
each of product form 

K t (y 7T1 5 Vn 2 ) ■ ■■Kt{y 7T2n-l ? IJ^n) (18) 
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for some permutation n, containing each of the variables (y,; : i < 2 n) exactly once. Hence 
each term, and therefore the entire Pfaffian, solves the desired equation dtu = Yl'f=i L yi v, 
when y £ V 2 

Note that \Kt(x,y)\ < 1 and hence the Pfaffian Pf {K^ 2 n \t,y)) is a uniformly bounded 
function on [0,oo) x Vi n . Uniqueness of solutions to the sequence (((ODE) 2n : n = 1, 2,...) 
now implies that 

E, ; [£ y (rit)} = u^ 2 n \t,y) = Pf (/\ (2n) (t,y)), for n > l,t > 0 and y G V 2n . 

This lemma states that this identity holds also for y £ Vin • However, by repeating the 
argument for the boundary conditions, for y £ V 2n any equalities in y\ < y 2 < ... y 2n can be 
removed, pair by pair, until 

u ( 2 n \t,y) = u^ 2 m \t,z), and Pf (K< 2 n \t, y)) = Pf {K^ 2 m \t,z)) 
for some m < n and z G V 2m , and hence equality also holds on the larger set. I 

Proof of Theorem [Tj The correlation functions E[ry(xi)... ry{x n )] can be recovered from 
the product spin expectations. Indeed 

D t ff yAv) = a v,z+i (v) ~ VyA 1 !) 

= <r y M (i-e)"™ - i) 

= -(l + 9)ri(z)<Ty iZ (rj) 

so that 

D t°yA 7 i) \ z=y = -(i + o)v(y)- ( 19 ) 

Thus out of a spin pair, we can reconstruct a single occupancy variable by first a discrete 
derivative, and then an evaluation. We now iterate this to get at multiple occupancy variables. 
Fix n > 1 and consider y = (aq, x \,..., x n , x n ) ^V 2n where 

X\ < X\ < x 2 < x 2 < X 3 < ... < x n < x n . 

The restriction that Xi < Xj+i allows us to apply the operators D, t,..., to both sides of 
the identity (from Lemma |h]) 

E, [S,(y t )]=Pf {K^\t, y)). 


The left hand side becomes 


(-!)"(! +0) n E„ 


n 

r]t{x i) • ■ -r]t(x n ) Y[ a x i ,x i ( r 1t) 

i —1 


After setting X{ = x* for all i we reach the correlation (—l) n (l + 0)"E^[y t (xi)... ?y(x n )]. 
Applying the operators D, t,..., D ,t to the Pfaffian on the right hand side preserves the 
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Pfaffian structure. Indeed applying to the single Pfaffian term (18) will change only a 
single factor in the product, namely 


Kt(yn2i-i 


Vl T 2 i) 


Df if ir 2i -i = 1, 

D ~2 K t {y n2i - 1 , y-K 2 i) if 7T2i = 1- 


The terms can then be summed into a new Pfaffian where the entries in the second row and 
column, which are the only entries containing the variable x±, are changed. Repeating this 
for the operators D t ,..., , and then setting Xi = x t for all i we still have a Pfaffian, and 

the hnal entries are in the form 

- 1 ( K t (y,z) D ~2 K t (y, z) \ 

1 + d V D+K t (y, z) D+D+K t (y, z) ) 


A final conjugation with a determinant one diagonal matrix, with entries that are alternating 
+i, will adjust the minus signs to give exactly by the kernel K stated in the theorem. I. 


Remarks. 

1. For random initial conditions the point process rjt will not i n general be a Pfaffian point 
process, although by conditioning on the initial condition, the intensities can always be 
written as the expectation of a Pfaffian. However, under some random initial conditions, ly 
does remain a Pfaffian point process. Indeed, examining the proof of Lemma [6] , one needs 
only that the expectation E [E y (r/o)], for y 6 V 2n . can be written as a Pfaffian Pf (<I yj) : 
i < j) for some $ : V 2 —> M. One then replaces the initial condition in the equation ([5]) for 
Kt(y , z) by &(y, z) and the rest of the argument goes through. An example is where the sites 
(rjo(x) :i£Z) are independent with r/o(x) a Bernoulli (p x ) variable, then 

E [£ 2 ,( 770 )] = Pf ( If (!-(! + 6 )p k ):i<j 

\k£{yi,yj) 



2. A slightly more combinatorial way of writing out the argument for the last part of the proof 
of Theorem [T] is as follows. Starting from (14) we may reconstruct the product intensities as 


E 




_i=1 


= E 


n 

_i =1 


(1 &Xi,Xi + l( 7 lt)) 

1 + 9 




E 


771=0 


y\<...<ym 


a yi,yi+^i r lt 


. 1=1 


( 20 ) 


Since the vector y( 2m ) = 
that 

n 


E 


Ylvt(xi) 
_i= 1 


( 2 / 1 , 2/1 + 1,..., 2/m, y m + 1) e v 2m we may apply Lemma 0 to see 

n 

= (l + 6)~ n Pf ^K< 2m \t,y (2m) )^ . 

m =0 Wl <-..<!Im 

e{xi,...,x n } 
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This sum may be recombined as the single Pfaffian 


E 


}~[vt( x i) 


_i= 1 


= (—!)”(! + 0)" n Pf (Kd 2n \t,y^) - J 2n ) 


where J 2n is the block diagonal matrix formed by n copies of J ). Indeed this is a special 
case of the general formula for Pf (A + B ) (see (31) in the appendix). This shows that rjt is 
a Pfaffian point process with the kernel K given, for y < z, by 


KG/,*) 


-1 

1 + 0 


Kt(y,z) K t (y,z + 1) \ 

K t (y+l,z) K t (y + l,z + l) ) ’ 


and Ki 2 (y, y) = j—y + 1) — 1), and other entries determined by the symmetry condi¬ 
tions. That both kernels K and K determine the same point process can be seen by a simple 
transformation: conjugation by an elementary matrix which subtracts the first row and col¬ 
umn from the second row and column, followed by a conjugation to adjust the positions of 
the minus signs as before, will transform the entries of K into those for K. 


3 Examples 


We prove the various continuum limits given as the four examples (A),(B),(C) and (D) 
listed in Theorem [2] The fact that the discrete kernels satisfy two dimensional discrete heat 
equations make all the examples easy to guess, and they are not that much harder to prove. 
We divide our proof into three steps: (1) convergence of the lattice particle systems to their 
continuum analogues for finite particle systems; (2) establishing the Pfaffian kernel for these 
finite continuum system by limits of the lattice Pfaffian kernels; (3) taking limits to obtain 
the continuum systems under maximal entrance laws and solving explicitly for the kernel. 
(In [22] we argued directly with the continuum Brownian models, and we don’t doubt that 
this would also be possible for these four examples.) 

We use the space Mlfp of locally finite point measures on M with the topology of vague 
convergence, which has simply checked compactness properties. Our point processes can then 
be thought of taking values in the measurable subset Mo C Mlfp of simple measures, that 
is where all atoms have mass one. 


Step 1. Finite particle continuum approximation. Working first with finite particle 
systems avoids complicated weak convergence arguments for infinite systems. We fix IV > 1 
and initial particle positions a\ < a 2 < ... < a n- In example (C) we assume ai > 0 and in 
example (D) we assume ai > 0. We choose finite lattice initial conditions 


Vo{x) 


1 if x = df- 1 for i = 1,..., N, 
0 otherwise, 


( 21 ) 
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where edf ' 1 —y m as e —> 0 for i = 1,..., N. We fix t > 0 and set x[ e \dx) = Ve~ 2 t{ x )^ex- 
The conclusion we want is that 

(X { t e) :te[0,T])^(X t :t€[0,T]) 

in the Skorokhod space D([ 0 , T],Mlfp)- One way to check this is via the continuous mapping 
principle (see PD, from the weak convergence of the scaled non-interacting simple random 
walk paths (Z*’ 1 ,..., Z*> N ) -> in D([0, T], W N ), to independent Brownian 

motions. For the purely annihilating system there is a deterministic map F : Z)([0, T], M^) —> 
D([0, T], Mlfp) that ’prunes’ the paths appropriately at collision times, satisfying 

(X[ e) :te[0,T}) = F(Zf\ ..., Z/ N : t E [0, T]). 

We will not detail the natural definition of F. but it is straightforward when the collision 
times of the N paths are all distinct, and and when two collision times are equal we may 
define F to be the zero function. The law of the limiting Brownian paths does not charge 
the discontinuities of F, and desired result follows. 

For mixed models the argument can be adjusted to allow for the random reaction that is 
required, and we give an informal description. From N initial particles there are at most 
N reactions. We can fix a vector of reactions r = (n, ...,r n) where r,; E {A,C + 
reflect the successive decisions: annihilate, coalesce onto the ’higher’ of the two colliding 
paths, coalesce onto the ’lower’ of the two colliding paths. For a fixed r there is again a 
deterministic map F r : D([ 0 , T], M. N ) —> D([0,T], Mlfp) that prunes the N paths according 
to the set of decisions (not all the decisions may be needed by time T). Then, for a bounded 
and continuous : D([0,T],Mlfp) ^ 

E[<5>{xl e) : t E [0, T])] = ^ 0 A{r \l - 9) N ~ A{r) E[^(F r (Z ^,..., Z e t ' N : t E [0,T]))] 

r 

where A(r) is the number of annihilations in r, and the result follows as for the coalescing 
case. For example ( D ) we apply the same to reflected paths, and for example (C) we include 
the times when a particle first hits zero as an extra collision, which reduces the number of 
particles by one but which requires no random reaction. 

Step 2. From lattice kernels to continuum kernels We record the limiting continuum 
Pfaffian kernels, for systems with finite initial conditions. 


Proposition 7 The continuum models with finite initial conditions r]Q = S ai are Pfaf¬ 
fian point processes at a fixed t > 0 with kernel in the form (Q) where (y, z ) is the unique 

solution to the following PDES: 

• (A),(B) Coalescing/annihilating Brownian motions. 

dt K t C) (y, z) = A k[ c) (y , z) for y < z, t > 0, 

< K$ c \y,y) = 1 for y E M, f > 0, (22) 

k l4 c \y,z) = (-6 /fory<z, 
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•(C) Coalescing/annihilating Brownian motions on [0, oo) killed at {0}. 


dtK[ c \y,z) 
Kl c \y,y ) 
diKf\0,z) 
, Kq\v,z) 


A K[ c) (y, z) for 0 < y < z, t > 0, 

1 for 0 < y, t > 0, 

0, for 0 < z, t > 0, 

(— Q/o[y,z) j or 0 < y < z. 


(23) 


• (D) Coalescing/annihilating reflected Brownian motions on [0,oo). 


d t K^ c \y,z) 

t K l c \y,y) 

hf\0,z) 
. K$\y,z) 

where (3>*(z) : z,t > 0) solves 


A K[ c) (y, z) for 0 < y < z, t > 0, 
1 for y > 0, t > 0, 

$t(z) for 0 < z,t > 0, 

f or 0 < y < z, 


f d t $t(z) = A $ t (z) for 0 < z, t > 0, 
1 ®t(0 ) = 1, fort> 0. 


(24) 


In each case we will show the convergence of the approximating lattice system — > X. 

The following lemma (whose proof is at the end of the appendix) is a natural approximation 
lemma for lattice kernels to continuum kernels. 


Lemma 8 For e > 0, let X^ e \dx) be random point measure on M whose atoms from a 
Pfaffian point process on eZ with kernel K^. Suppose that 

sup ||e _ 1 K^)|| 0O = sup sup e -1 |K^(y, z)\ < oo (25) 

e>0 e>0 y,z,i,j 


and 


lime 1 K^ ) (j/ e , z e ) = K \f(y, z), for i,j e {1, 2}, 
e4-0 J J 

when ( y e , z e ) -» (y, z) with y < z, or when y e = z e ^ z = y, 


(26) 


for some continuum kernel K (c) ■ r2 R 2x2_ Then X (e) X 

in distribution as e | 0, on 

the space AAlfp and the limit X is simple, and is a Pfaffian point process with kernel K^ c h 


Note that in our examples the limiting kernel K^(x, y) will be discontinuous at x = y. 


(A),(B) . This is the case p x = q x = 1 for all x E Z. The corresponding one-particle generator 
is then the discrete Laplacian A = D + + D~~. Lemma [5] shows that the entries of the Pfaffian 
kernel for X^ are given for y < z in eZ by 


K U(y,z) 


l 

1+0 


K e - 2 t (e ^,6 ^ -D%K e - 2 t (e 1 y,e l z) \ 
—Df K e - 2 t (e~ 1 y, e~ x z) Df D/K e -i t (e~ l y, e~ 1 z) ) ’ 
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and by K^(y,y) = ^ D^K e - 2 t (e } y,e 1 y). 


We rescale the scalar kernel by defining 


Kt e \y, z ) = K e - 2 t (e 1 y, e l z) for y, z <E eZ. 


For / :eZ 


we set 


D&>+f( X ) = f( x + e )~f( x \ A^f(x) = f( x + e ) + f( x - e )~ 2 f( x \ 


Then 


K^\y,z) = 


Kf\y,z) - eDi e) ’ + hf\y,z ) 

-eZ^’+A^y, z) e^D[ e) ’ + D^ + hf\y, z ) 


and Kf 6 \y, z) solves, for y,z E eZ, 


^AT t (e) (y, z) = (A^ e) + Ai e) ^) A t (e) (y, z) for y < z, t > 0 , 

‘ K ( t e \y,y) = 1, for all y, t > 0, 

I< { 0 e \y,z) = (- 0 )’ 7 o[ £ - 1 2 /, 6 -A) for y < z. 


(27) 


(28) 


By conjugation with a diagonal matrix with diagonal entries e 1 / 2 and e x / 2 an alternative 
kernel for X^ e \ in the right form for Lemma pj is 


K {e \y,z) = 


—D^' + K^ e \y,z) 

T+ 0 \ -D®' + Kj: e \y,z) D[ e) ’+D^ +t Ki e \y,z) 


K t\y,z) 


Checking the hypotheses (25) and (26) amounts to checking that the lattice approximations 
to the two dimensional continuum PDE (22) converge uniformly, at a fixed t > 0, along with 
their first and second derivatives. The required estimates are quite standard and we omit the 
proof here, and also for the examples C,D below. Some details (however for different initial 
conditions) are contained in the thesis 


(C). We take 


1 for x > 1 , 
q x = { 2 for x = 0 , 
0 for x < 1 , 


Px = 


1 for x > 1 , 
0 for x < 1 . 


(29) 


No particle ever visits {..., —3, — 2 }, and particles that reach {—1} never escape. So we re¬ 
strict attention to the point process (rjt(x) : x = 0, 1 ,...) and define X^ ( dx ) = X)z>o r le- 2 t( x )^ex 
as a measure on eZ + = {0, e, 2e,..A single particle acts as a simple random walk on 
{0,1,...} with a certain rate of being killed whenever it is at zero. Under diffusive rescaling 
this process becomes a Brownian motion that is instantly killed at the origin (which allows 
the treatment for interacting particles to go through as in step 1 .) 

The reason for choosing qo = 2 is so that the corresponding one particle generator L defined 
in Q is the generator for a reflected random walk on { 0 , 1 ,...}, which jumps x —> x ± 1 at 
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rate 1 for x > 1, and also jumps 0 —> 1 at rate 2. This can be realised as the absolute value of 
a simple random walk on Z, and this helped us in some of the estimates showing the lattice 
PDE converged to the continuum PDE. 


Then is Pfaffian with kernel of the form |27| where Lemmajsjshows that K^ e \y, z ) solves 
(28) for e < y < z. As expected, the Neumann boundary condition emerges in the limiting 
continuum PDE (23) when e —> 0. 


(D). To obtain reflected random walks on Z + we take 


f 1 for x > 1, _ f 1 for x > 1, 

\ 0 for x < 1, j0,r | 0 for x < 1. 


We may restrict attention to the process ( rjt(x ) : x > 0). Note that the corresponding one 
particle generator L defined in Q is the generator for simple random walk on {0,1,...} 
absorbed at 0. Thus is Pfaffian with kernel of the form (]2T|) where (y, z) solves 


for e < y < z. The boundary condition can be found by examining the expectation 
K t ( 0, z) = [o-q , z {r}t)\- Applying the generator we find that 


d t K t (0, z) = L z K t ( 0, z) for z > 1, 


with the boundary condition Kt( 0,0) = 1. Under scaling the limiting boundary condition 
becomes &t( z ) as stated in example D of Theorem [2j 


Step 3. Maximal entrance laws. Infinite systems of coalescing particles are easy to build 
by adding one particle at a time. Infinite systems of annihilating particles perhaps require 
a bit more care. In [22] a pathwise construction is bypassed since only fixed time properties 
are studied, and instead this paper uses a Feller transition kernel pt(p, du) on the measurable 
subset Mo of simple measures within the space Mlfp • Our mixed models also have Feller 
transition kernels, and we can construct these by exploiting the Pffafian structure. Suppose 
—> r] E Mo, where 7/00 have finitely many particles. Then 

(_0^ N) [y,z) ^ (_0)ijb/,*) for a i mos t a ll z ). 


Indeed this holds for all y,z that are not the position of atoms in ?/. These functions are 
the initial conditions for the pdes (22), (23), (24) that determine the Pfaffian kernels. The 
solutions to these pdes then converge, at a fixed t > 0, together with their first and second 
derivatives. Thus the associated Pfaffian kernels converge (in a bounded pointwise manner) 
and Lemma [ 9 ] in the appendix states this is sufficient for the associated point processes 7/1^ 
to converge in law to a limit. This defines a kernel ptirpdv) for all r) E Mo- To check the 
semigroup property it is sufficient to have the Feller property, and pass to the limit from the 
semigroup property for finite systems. The Feller property follows however from the same 
argument: the convergence —> 7/ E Mo implies, via convergence of the Pfaffian kernels, 

the weak converge of the laws pt(j]^ N \du) —> Ptiry dv). 


To construct the maximal entrance laws we pick any sequence rj' N ' so that 


(_ 0 )U Ar) [yU ^0 as A -> 00 . 


(30) 
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When 0 G [0,1) it is easy to ask for pointwise convergence, but when 9 = 1 this is impossible. 
Instead we ask that the convergence holds in distribution (either on M or on [0, oo) for the 
killed and reflected models). This is still sufficient to imply that the associated PDEs, and 
its derivatives, converges pointwise at a fixed time t > 0 , as can be verified by writing out the 
solution in terms of the initial conditions and the associated Green’s functions. The limiting 
laws we denote as pt( oo, du), and again the Feller property allows one to check that they act 
as entrance laws for the Markov family. 


The name maximal is natural for coalescing systems, while for purely annihilating systems 
many increasingly dense initial conditions will not converge to pt{oo,dv) - a sequence of 
closely positioned pairs will annihilate quickly and leave empty regions. The Pfaffian structure 


makes it clear that the convergence in distribution in (30) is exactly the right condition to 


describe the domain of attraction for the maximal entrance law. The convergence (30) holds 


in particular for lattice initial conditions as the spacing decreases to zero, or by choosing 
Poisson initial conditions of increasing intensity. 


Finally under the maximal entrance laws the associated pdes (22), (23), (24) have zero initial 


conditions, and the explicit solutions listed in Theorem [2] are straightforward to derive. 


4 Appendix: Pfaffian Facts 


We collect here some basic facts about Pfaffians. All the following results are contained, for 
example, in Stembridge M- 

Throughout this section A = l is a 2n x 2 n anti-symmetric matrix (say with complex 

entries). The Pfaffian can be defined by 

1 n 

Pf( ^ )= 2 ^d ^ S gnWn°— 

7ren 2 „ i =i 

where IFn is the set of permutations on {1, 2,..., 2?z} and sgn(7r) is the sign of the permuta¬ 
tion 7 T. 

Pfaffians are related to determinants and determinant properties often have Pfaffian ana¬ 
logues. Many of these results are consequences of the following conjugation formula 

Pf ( CAC t ) = det(C) Pf (A) for all 2 n X 2 n matrices C, 

and the identity det(A) = (Pf (A)) 2 . 

One may decompose a Pfaffian along any row (or column) in terms of sub-matrices, analogous 
to the Laplace expansion of a determinant: for any i E {1,..., 2n} 

2 n 

Pf(A)= XI (_l)i+i+i+x{i<d a . jPf 

3= 


22 








where is the (2 n — 2) x (2 n — 2) sub-matrix formed by removing the i-th and j-th rows 
and columns from A. 

For two anti-symmetric matrices A, B, we have 

Pf (A + B) = ^(-i)l^l/ 2 (- 1 )^) P f (A\ v )Pf (B\ uc ), 
u 

where the sum is over subsets U C {1,..., 2n} with \U\/2 e {1,..., n}, s(U) = Ylj^u■) ( w ith 
s(0) = 0), and A\u means the matrix A restricted to the rows and columns of U. 

We use the special case of this when —B = J 2 n , the canonical symplectic matrix consisting of 

n blocks ^ ^ ^ ^ down the diagonal. Then setting U x = {2x\ — 1, 2x \,..., 2x m — 1, 2x m } for 

x € {1,2,... , n} m with x\ < X 2 < ■ ■ ■ < x m , we have Pf ( J 2 n\U x ) = 1 while Pf ( J 2 n\U ) = 0 
for all (non-empty) U C {l,2,...,2n} not of this form. Note that (— = (— l) m . In 

this case, the formula reduces to 

n 

Pf (A - hn) = XI E(- 1 >”"”' Pf f- 4 !")’ ( 31 > 

m =0 x 

where the sum over x is precisely the sum over x of the above form. 


Proof of the Pfaffian identity ( |17[ ). Let A be the antisymmetric matrix with entries 
A{j = CLi/aj. By conjugating with a suitable elementary matrix we may subtract a multiple 
01/02 of the second row and column form the first row and column. This produces a new 
matrix A with Pf (A) = Pf (A) but where the top row of A is now (0, ai/ 02 , 0,0,...). By a 
Laplace expansion of this top row we find 

Pf (A) = Pf (A) = —Pf = —Pf (^4 (1 ’ 2) ), 

Ci 2 CL 2 

and by induction over n we find Pf (^4) = (0103 ... 02 n-i )/(«204 • • • « 2 n)- I 


Proof of Lemma [8| For a < b, the first moments 


E 


* (e) (M) 


K[^(x,x) 

irGeZn [a, b] 


are uniformly bounded by (25), and this implies the tightness of ( X^ : e > 0) as elements of 
. Mlfp ■ The higher factorial moments (writing [z\k = z(z — 1)... (z — k + 1)) such as 


E 


[V ( ‘>([o, b])h 


= Xi Pf (K ( ^(xi,Xj) :i,j < k) 

iri,...,a:/ c GeZn [a,6] 

< ||e- 1 K^||^(( 6 -o) + e- 1 ) fc 
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are also uniformly bounded in e E (0,1]. By Fatou’s lemma, the moments of any limit points 
satisfy E [[V p )([a, &])]&] < C^fb — a) k for all k > 1. This implies that any limit point is a 
simple point process: indeed 


P [V({x}) > 2, for some x E [—M,M]\ 


< 


< 


< 


Mm 

J2 P[X([k/m,(k + l)/m)])>2] 

k=—Mm 

Mm 

£ E[[X([t/m,(t + l)/m)])] 2 ] 

k=—Mm 

Mm 

C 2 22 m ~ 2 —> 0 as m —> oo. 

k=—Mm 


Then one can pass to the limit, for finite disjoint intervals A\,..., A r 


E 




lk=\ 


22 pf (K {c \xi,Xj) :i,j < k) 


Xi^AiDeZ: i<k 


I Ai X...X Art 


Pf (K (c ) (. Xi,Xj ) : i,j < k)dx\... dx^ 


= E 


II X(A k ) 


lk=\ 


for any limit point X. Indeed the convergence follows from the assumptions (25) and (26); 
since E[X({x})] = 0 for limit points the discontinuities of the function fi — > n{A{)... n(A m ) 
are not charged; the finite higher moments give the uniform integrability that justify the final 
equality. This shows that a limit point is a Pfaffian point process with kernel K^-*. Finally 
it is well known (see the remark after this proof) that the fact that the kernel is locally 
bounded is sufficient to determine the law of the associated Pfaffian point process in Mlfp, 
which implies that limit points are unique and thereby the convergence of X^ e \ I 


Very similar arguments establish the following kernel convergence lemma for continuous ker¬ 
nels. 


Lemma 9 Suppose (X^ N \dx) : N > 1) are random point measure on M whose atoms form 
Pfaffian point processes on M with kernel Suppose that 

sup IIK^Iloo = sup sup \K.\P(y,z)\ < oo 

N N y,z,i,j 

and 

lim K \f\y,z) = K ij(y,z), fori,j E {1,2}, andy,z eR 

N^oo J 

for some limiting kernel K : M 2 —> R 2x2 . Then X^ —> X in distribution as N —^ oo, on the 
space Mlfp cmd the limit X is simple, and is a Pfaffian point process with kernel K. 
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Remark. To see that a locally bounded kernel K is sufficient to determine the law of 
associated point process X one may argue as follows: (i) all moments of X(f), where 
/ = e;=i cii(Ai) > 0 and Ai are finite intervals, are given in terms of the kernel; (ii) 
by Hadamard’s inequality 

|Pf (K (xi,Xj) : i,j <k) | = | det(K(xj, xj) : i,j < k) | 1/2 < ||K||^(2A:) fc ; 

(iii) the factorial moments E[[X([a, 6])]*] are bounded by C(a, 6 , K) fc /c fc ; (iv) the moments 
E[|X([a, 6])| fc ] are also bounded by C'(a,b,K.) k k k ; (v) the moment problem for X(f) is well 
posed. 
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